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ABSTRACT 


Adaptive,  frequency  domain,  array  signal  detection  is  considered  when  there  exist  data  blocks 
containing  interference  alone  and  signal  plus  interference  that  are  disjoint  in  time  and  span  several 
frequencies.  The  data  block  solely  containing  interference  is  used  to  form  estimates  of  the 
interference  covariance  structure  at  each  frequency.  These  estimates  are  substituted  into 
generalized  likelihood  ratios  to  form  constant  false  alarm  rate  detectors  for  deterministic  and 
Gaussian  signals  that  are  limited  in  time  to  the  signal  plus  interference  data  block.  Central  limit 
theorem  based  normal  approximations  are  used  to  determine  thresholds  and  the  signal-to- 
interference  ratio  (SIR)  required  to  achieve  specified  false  alarm  and  detection  probabilities.  The 
SIR  required  for  the  proposed  detectors  is  shown  to  be  the  SIR  required  for  the  ideal  conventional 
detector  with  interference  covariance  estimation  losses  and  intra-block  correlation  losses.  The  gain 
in  adaptive  processing  (i:e.,  the  reduction  in  signal  strength  required  to  achieve  specified  false 
alarm  and  detection  probabilities)  is  seen  to  be  the  ratio  of  the  array  gain  improvement  (AG/)  to  the 
change  in  the  estimation  and  correlation  losses  resulting  from  the  shift  from  conventional 
beamforming  (dimension  equals  one)  to  higher  dimension  processing.  Thus,  the  AGI  must 
exceed  the  change  in  the  estimation  and  correlation  losses  before  adaptive  detection  becomes 
attractive. 
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A  MAXIMUM  LIKELIHOOD  BASED  ADAPTIVE  DETECTOR 

1  INTRODUCTION 

Adaptive  detection  of  signals  received  by  an  array  of  sensors  provides  improved  performance  over 
contemporary  methods  involving  conventional  beamforming  prior  to  a  detection  algorithm  in 
interference  environments  that  yield  high  array  gain  improvement  ( AGI ).  Several  issues  plague 
adaptive  detection,  the  foremost  involves  the  estimation  of  unknown  covariance  matrices. 
Imperfect  knowledge  of  covariance  matrices  both  degrades  performance  and  complicates  the 
detector  derivation  and  analysis.  Detectors  are  often  based  on  optimal  techniques  drawn  from  the 
Neyman-Pearson  theory  of  hypothesis  testing.  Optimal  in  this  case  implies  maximizing  the 
probability  of  detecting  a  signal  given  a  fixed  probability  of  false  alarm.  When  a  complete 
statistical  description  of  the  signal  and  interference  is  available,  these  techniques  provide  the  best 
detector  structure  in  the  form  of  a  likelihood  ratio  test  (LRT).  When  only  partial  information  about 
the  statistics  is  available,  optimal  tests  rarely  exist  and  the  suboptimal,  although  usually  very  good, 
generalized  likelihood  ratio  tests  (GLRT's)  are  used.  Ideally,  a  detector  incorporates  known  signal 
information  such  as  time  and  frequency  characteristics  into  the  design  of  the  algorithm.  This 
matching  of  the  detector  structure  to  the  signal  unfortunately  means  that  a  loss  is  incurred  if  the 
signal  does  not  have  the  assumed  structure,  perhaps  due  to  time  or  frequency  spreading  of  a 
transmitted  signal,  propagation  effects  of  the  medium,  or  simply  a  lack  of  knowledge  of  signal 
characteristics.  A  detector  that  is  robust  against  not  only  signal  mismatch,  but  also  against  the 
general  time-frequency  characteristics  of  the  signal,  is  desired. 

It  is  commonly  assumed  that  the  interference  hindering  detection  is  stationary  over  a  length  of  time 
greater  than  that  of  the  signal.  In  frequency  domain  array  processing  (i.e.,  discrete  Fourier 
transform  (DFT)  of  the  array  time  data),  the  signal  may  be  assumed  to  lie  in  a  block  of  data  in  time 
(consecutive  DFT  outputs,  also  called  time  snapshots)  and  frequency,  denoted  the  primary  data 
block  (PDB).  Under  the  above  stationarity  assumption,  the  interferences  present  in  the  PDB  are 
also  present  in  data  in  adjacent,  disjoint  time  snapshots  at  the  same  frequencies.  This  data  is 
grouped  and  denoted  the  secondary  data  block  (SDB).  The  primary-secondary  data  block  structure 
is  shifted  in  time,  continually  testing  for  the  presence  of  a  signal  in  the  PDB  subject  to 
normalization  by  the  interferences  present  in  both  data  blocks.  Wolcin  [1],  with  extensions  by 
Nuttall  [2],  has  considered  such  a  PDB  with  exponential  averaging  of  the  interference  covariance 
structure  over  a  semi-infinite  SDB  consisting  of  all  data  previous  to  the  PDB.  Specific  subsets  of 
these  general  block  assumptions  are  found  in  the  cardinal  work  of  Reed,  Mallett  and  Brennan  [3], 
as  well  as  more  recently  in  several  papers,  notably  by  Kelly  [4],  Robey  et  al.  [5],  Robey  [6],  and 
Chen  and  Reed  [7].  These  papers  [3-7]  consider  the  single  frequency  case  where  the  PDB  consists 
of  a  single  array  vector  and  the  SDB  is  a  set  of  vectors  with  the  same  interference  covariance 
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structure.  Robey  [6]  does  consider  the  case  of  a  multiple  vector  PDB,  however,  the  detector 
analysis  is  for  the  single  vector  case.  These  papers  are  also  based  on  a  signal  with  deterministic 
form  and  an  unknown  complex  amplitude. 

This  memorandum  provides  a  brief  background  on  array  signal  detection,  proposing  the  use  of 
matrix  preprocessing  to  reduce  the  adaptive  dimension  of  the  detector  and  discussing  common 
statistical  models  of  the  array  data.  Detectors  will  be  developed  and  analyzed  for  the  deterministic 
signals  considered  in  the  above  references  [3-7]  as  well  as  a  Gaussian  signal  with  unknown 
power.  The  detector  development  is  based  on  a  generalization  to  the  likelihood  ratio  approach 
reported  by  Robey  [5,  6].  Following  Wolcin  [1],  the  detector  for  the  Gaussian  signal  is  designed 
to  account  for  varying  signal  time-frequency  characteristics  by  only  considering  time-frequency 
cells  in  the  PDB  that  may  contain  a  signal  exceeding  a  particular  strength.  This  effectively 
estimates  the  location  of  the  signal  in  the  PDB  time-frequency  cells. 

The  detector  performance  is  evaluated  by  the  signal-to-interference  ratio  (SIR)  required  to  achieve 
desired  false  alarm  and  detection  probabilities.  When  there  are  enough  cells  in  the  PDB,  analysis 
may  be  performed  by  a  normal  random  variable  approximation  using  the  central  limit  theorem 
(CLT).  The  SIR  is  decomposed  into  the  required  SIR  for  the  ideal  non-adaptive  detector  (i.e., 
conventional  beamforming  prior  to  detection  when  the  conventional  beam  output  power  is  known) 
as  reported  by  Wolcin  [1]  or  Nuttall  [2],  with  losses  due  to  the  estimation  of  the  interference 
statistics  and  losses  due  to  intra-block  correlations  caused  by  the  estimated  interference  statistics. 
The  gain  in  adaptive  processing  (i.e.,  the  reduction  in  the  signal  strength  required  to  achieve 
specified  false  alarm  and  detection  probabilities)  is  seen  to  be  the  ratio  of  the  AGI  to  the  change  in 
the  estimation  and  correlation  losses  from  univariate  to  multivariate  data. 

Applications  of  the  proposed  detectors  include  any  limited  time  duration  signal  in  an  interference 
background  with  moderately  constrained  stationarity  conditions.  The  degree  to  which  adaptive 
processing  is  applicable  depends  on  the  severity  of  the  interference  background.  The  AGI,  which, 
for  instance,  is  large  for  directions  near  to  strong  plane  wave  (isolated  in  arrival  angle) 
interferences,  must  exceed  the  change  in  the  estimation  and  correlation  losses  before  adaptive 
processing  becomes  attractive. 


2  ARRAY  SIGNAL  DETECTION  BACKGROUND 

The  conventional  detection  of  a  signal  arriving  at  an  array  of  sensors  from  a  particular  direction 
may  be  viewed  as  two  operations.  First,  the  multivariate  array  data  is  processed  to  form  a 
univariate  test  statistic.  This  is  often  accomplished  by  conventional  beamforming  and  is  discussed 
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in  section  2.1.  Second,  a  function,  called  the  detector  function  or  detector  non-linearity,  is  applied 
to  the  univariate  test  statistic.  The  detector  function  is  formed  based  on  the  statistical  characteristics 
of  the  test  statistic  which  are  dependent  on  the  statistical  assumptions  on  the  array  data.  The  output 
of  the  detector  function  is  then  compared  to  a  threshold  to  decide  if  a  signal  is  present  or  not. 

Adaptive  array  detectors  form  a  detector  function  based  on  the  statistical  characteristics  of 
multivariate  array  data  rather  than  the  univariate  processed  data.  Unfortunately,  this  method  is 
fraught  with  analysis  and  implementation  difficulties.  The  estimation  and  inversion  of  the  array 
data  covariance  matrix  is  typically  required.  As  the  number  of  array  sensors  increases,  the  matrix 
inversion  becomes  numerically  inhibitive  and  the  matrix  estimation  becomes  statistically  hazardous. 
A  method  of  dealing  with  large  arrays  involves  preprocessing  the  array  data  to  a  smaller  dimension 
and  then  applying  an  adaptive  detection  algorithm.  A  very  general  matrix  preprocessor  is  described 
in  section  2.1,  followed  by  a  specific  method,  beam  space  preprocessing,  that  is  applicable  to 
adaptive  detection.  Common  statistical  models  for  the  array  data  and  the  effects  of  the  matrix 
preprocessing  on  these  models  are  discussed  in  section  2.2. 

2.1  Pre-Detector  Signal  Processing 

The  signal  processing  of  an  array  of  sensors  in  sonar  or  radar  detection  applications,  as  seen  in 
figure  2- 1,  typically  includes  simultaneous  time  sampling  of  the  continuous  waveform  at  each  of  P 
sensors  followed  by  transformation  to  the  Fourier  domain  by  a  discrete  Fourier  transform  (DFT). 
The  frequency  domain  array  data  is  beamformed  to  multiple  angles  with  a  signal  detection 
algorithm  at  each  beam  output. 

Beamforming  is  traditionally  performed  to  increase  the  signal-to-noise  ratio  (SNR)  prior  to  the 
detector  by  spatially  filtering  noise  and  interferences  from  other  angles  of  arrival.  This  may  be 
accomplished  by  using  conventional,  non-adaptive  beamforming,  or  by  data  dependent,  adaptive 
methods. 


DFT 

-  \ 

Beamformer 

— 

Detector 

—  / 

frequency  fm 

= 

angle  0( 

Figure  2-1  Typical  signal  processing  of  sensor  arrays  for  detection. 
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A  more  general  configuration  of  this  system  replaces  the  beamforming  block  with  an  angle 
dependent  matrix  preprocessor  followed  by  a  detector  as  seen  in  figure  2-2.  If  the  array  data  for 
the  m,h  DFT  bin  and  the  Ith  time  snapshot  is  arranged  in  the  P-by- 1  vector, 


Xm/  —  ml  *2 ml  ' '  '  */wl 


(1) 


and  operated  on  by  the  P-by- A  preprocessing  matrix,  AmS,  the  resulting  TV-by-1  preprocessed 
data  vector  is 


(2) 


where  the  superscript  T  is  the  transpose  operation  and  the  superscript  H  is  the  conjugate  transpose 
operation.  Vectors  are  denoted  by  lower  case  bold  letters,  and  matrices  are  denoted  by  upper  case 
bold  letters.  The  "  notation  found  on  vector  variables  indicates  that  the  variable  is  at  the  array 
sensor  level. 
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Figure  2-2  Preprocessing  and  adaptive  detector. 

Preprocessing  the  array  sensor  data  provides  V-dimensional  data  to  the  detector,  where  N  is  now 
dependent  on  the  type  of  preprocessor  employed.  Conventional  beamforming  is  represented  by 
setting  the  matrix  preprocessor  to  the  array  steering  vector  for  the  desired  angle  of  arrival. 
Similarly,  element  space  adaptive  processing  is  represented  by  setting  the  matrix  preprocessor  to 
the  P  dimensional  identity  matrix,  thus  passing  all  of  the  sensor  data  to  the  detector.  The  matrix 
preprocessing  structure  also  allows  reduction  of  the  adaptive  dimension  by  transformations  such 
that  1  <  N  <  P,  as  in  the  subarray  and  beamspace  preprocessors  discussed  in  [8]. 

Beamspace  preprocessors,  as  seen  in  [8],  are  of  particular  interest  because  the  adaptive  dimension 
may  be  greatly  reduced  with  minimal  performance  loss. .  For  example,  the  AG1  for  an  array  with 
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P  =  64  sensors  using  N  =  5  beams  in  a  beamspace  preprocessor  is  typically  within  0.5  dB  of  the 
fully  adaptive  beamformer  AGI.  The  preprocessing  matrix  has  the  form 

A.»=[d.,(e)  ...  a_  <«].  (3) 

where  dmi(#)  represents  an  array  steering  vector  at  the  m'h  DFT  frequency  with  a  direction 
denoted  by  the  index,  i.  The  N  beams  that  are  chosen  are  those  equally  spaced  in  wavenumber 
and  centered  about  the  desired  angle,  6. 

When  adaptive  beamforming  (and  also  detection)  is  employed,  reduction  of  the  adaptive  dimension 
is  crucial  due  to  the  numerical  intensity  of  the  implementation  and  the  statistical  hazards  involved  in 
the  estimation  of  large  covariance  maU'ices.  The  reduced  dimension  adaptive  techniques  typically 
require  estimation  and  inversion  of  an  N- by-  N  covariance  matrix  at  each  DFT  bin  output.  From  a 
numerical  computation  and  a  statistical  estimation  standpoint  this  becomes  difficult  as  N  increases. 
Preprocessing  the  array  data  vectors  to  provide  small  N  reduces  numerical  computations  often 
with  near  equivalent  average  performance  as  shown  by  Owsley  and  Abraham  [8]  and  also  provides 
improved  statistical  performance  as  shown  by  Burgess  and  VanVeen  [9]. 

2.2  Statistical  Modeling 

The  array  data  vector  from  the  m,h  DFT  frequency  bin  and  the  l,h  time  snapshot,  xm/,  is  commonly 

modeled  as  a  complex,  P  variate,  Gaussian  random  vector.  The  following  notation  will  be  used 
throughout  this  paper  to  represent  a  P  variate,  complex  Gaussian  random  vector  with  mean,  p., 
and  covariance  matrix,  X : 

x  ~  CAV(p,Z). 

When  no  signal  is  present,  denoted  by  the  hypothesis  H ,  the  array  data  vector  is  distributed  as 

xjtf  ~av,(0,Q.),  (4) 

where  Qm  is  the  P  -by-  P  covariance  matrix  representing  the  correlations  between  sensors  at  the 
m,h  DFT  frequency  bin  due  to  interference  and  noise. 

There  are  two  common  signal  models,  the  stochastic  Gaussian  signal  with  unknown  power  (SG) 
and  the  deterministic  signal  with  unknown  complex  amplitude  (SD).  Both  consider  a  plane  wave 
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arrival  with  either  stochastic  or  deterministic  complex  amplitude.  For  the  SG  case,  the  data  vector 
is  distributed  as 


ijK~CNP(Q,smll„d!,6  +  Qm)  (5) 

where  the  hypothesis  K  represents  signal  present,  dm8  represents  the  steering  vector  for  the  sensor 
level  data  and  smi  is  the  power  of  the  signal  in  the  m'h  DFT  bin  and  llh  time  snapshot.  The  SD 

case,  common  in  radar  array  processing,  considers  the  time  domain  signal  to  be  a  sinusoid  with 
known  frequency  and  unknown  complex  amplitude,  am, ,  yielding 

-  CNp(ccnildmg,Qm).  (6) 


For  comparison  purposes,  it  will  be  assumed  that  the  powers  in  the  deterministic  and  Gaussian 
signals  are  the  same,  i.e.. 


The  data  vectors  are  assumed  to  be  independent,  because  of  the  Gaussian  data  assumption,  for 
different  frequencies  and  different  time  snapshots.  Mathematically,  this  is  described  by 

Cov[xm,x%  ]  =  SmnS,jCov[xmli"  ] ,  (7) 

where  Sij  t  the  Kronecker  delta  function,  is  non-zero  only  when  i  =  j  where  it  has  unit  value.  The 

justification  for  this  is  based  on  [10]  where,  for  the  stationary  continuous  case,  Hodgkiss  and 
Nolte  have  shown  that  if  the  observation  time  of  the  signal  is  long  enough,  different  Fourier 
coefficients  become  uncorrelated  and,  thus,  are  independent  due  to  the  Gaussian  data  assumption. 
It  is  also  easily  shown  that  similar  results  hold  for  the  discrete  case;  that  is,  if  the  time  series  data  is 
stationary,  different  DFT  bins  are  approximately  uncorrelated  for  a  large  enough  DFT  size  and  thus 
independent,  assuming  non-overlapping  DFT's.  Note  that  since  our  observation  intervals  are 
short,  this  places  a  minimum  size  constraint  on  the  length  of  the  DFT. 


When  the  matrix  preprocessing,  described  in  section  2.1,  is  incorporated  into  the  statistical 
description  of  the  data,  the  SG  and  SD  signals  are  respectively  distributed  as 

X JK  -  +  Q„s )  (8) 
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and 

xm/l^  ~  CNw(ttm/dme,Qme ) 


(9) 


where  the  preprocessed  steering  vectors  and  interference  covariance  matrices  are 


and 


d/«e  —  Amed  mB 


Qn,e  =  AlQmA 


Due  to  the  generality  of  the  preprocessed  data  structure,  the  detection  problem  will  be  developed 
assuming  that  preprocessing  has  already  been  performed  and,  for  notational  convenience,  suppress 
the  dependence  on  the  angle  6 . 


3  MAXIMUM  LIKELIHOOD  BASED  ADAPTIVE  DETECTORS 

As  previously  mentioned,  the  array  data  is  segmented  in  time  into  the  PDB,  containing  interference 
and  potentially  signal,  and  the  SDB,  containing  interference  only.  This  structure,  the  statistical 
representation  of  the  primary  and  secondary  data,  and  its  potential  applications  in  the  field  of  sonar 
are  discussed  in  section  3.1.  Then,  following  Robey's  [5]  derivation  for  the  deterministic  signal 
case,  a  detector  for  the  Gaussian  signal  case  is  derived  in  section  3.2,  and  it  is  shown  that  each 
preprocessed  array  data  vector  in  the  PDB  is  used  to  form  an  adaptive  beam  output.  The  statistical 
distribution  of  the  adaptive  beam  output  is  then  described  in  section  3.3  for  Gaussian  and 
deterministic  signals. 

3.1  Array  Data  Structure 

The  primary-secondary  data  structure,  commonly  found  in  radar  systems  analysis,  is  composed  of 
two  sets  of  independent  data  blocks.  Each  data  block  is  composed  of  a  set  of  DFT  bin  outputs  for 
multiple  frequencies  and  time  snapshots.  The  primary  data  block  consists  of  data  that  is 
hypothesized  to  include  signal  energy.  The  secondary  data  block  is  independent  of  the  PDB  and 
must  have  the  same  interference  statistics  as  that  found  in  the  primary  data  block  at  each  frequency. 
The  SDB  is  used  to  form  estimates  of  the  frequency  dependent  interference  covariance  matrix. 
This  estimate  is  crucial  to  the  constant  false  alarm  rate  (CFAR)  property  of  the  proposed  detector. 
Figure  3-1  shows  a  typical  contiguous  primary  and  secondary  data  block  set  with  M  DFT 
frequencies,  L  time  snapshots  in  the  PDB,  and  J  time  snapshots  in  the  SDB. 

Each  time  snapshot-DFT  frequency  bin  point  in  the  data  blocks  represents  a  length  N  preprocessed 
data  vector.  Those  in  the  SDB  are  assumed  to  contain  no  signal  components,  where  those  in  the 
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PDB  may.  Statistically,  the  data  vectors  are  described  by  equations  (8)  and  (9),  respectively,  for 
the  SG  and  SD  cases.  The  no  signal  requirement  of  data  vectors  in  the  SDB  implies  that 
sml  =  0  or  aml  =  0.  This  yields  the  statistical  representation 


y «  -  cw„(o.q.)  (10) 

for  both  signal  types,  where  ymj  is  the  secondary  data  vector  from  the  m,h  DFT  frequency  and  the 
j'h  time  snapshot  (where  j  =  )  and  Qm  is  the  interference  covariance  matrix. 
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Figure  3-1  Primary  and  Secondary  Data  Block  representation. 

The  SDB  is  not  required  to  be  a  contiguous  set  of  time  snapshots.  Windows  leading  and  lagging 
the  PDB  may  be  used  to  estimate  the  interference  covariance  matrix,  as  in  cell  averaging  CFAR 
systems,  as  long  as  steps  are  taken  to  insure  that  these  windows  are  not  contaminated  by  signal 
presence  and  have  the  same  interference  covariance  structure  as  the  PDB. 

The  most  restricting  requirement  of  this  data  structure  lies  in  the  required  stationarity  of  the 
interference  over  both  the  primary  and  secondary  data  blocks.  The  detector  performance  is  nearly 
optimal  when  the  SDB  contains  a  large  number  of  time  snapshots;  this,  however,  conflicts  with 
typical  stationarity  assumptions  for  sonar  interferences.  The  requirement  that  the  SDB  not  contain 
any  energy  due  to  the  signal  limits  the  application  of  this  data  structure  to  scenarios  that  have  a 
definitive  signal  start  or  stop  time. 
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The  PDB  must  be  sized  based  on  the  type  of  signal  that  is  to  be  detected.  Ideally,  if  the  time- 
frequency  characteristics  of  the  signal  were  known  exactly,  only  the  cells  that  contained  signal 
energy  would  be  included  in  the  PDB.  However,  since  this  information  is  rarely  known  and  a 
detector  that  is  robust  against  unknown  signal  lime-frequency  characteristics  is  desired,  a  block  in 
time  and  frequency  is  utilized. 

These  restrictions  on  the  PDB  and  SDB  indicate  that  this  data  structure  may  be  applied  to  limited 
time  duration  signals  in  a  background  that  is  stationary  over  some  period  of  time  greater  than  the 
length  of  the  signal.  Included  in  this  category  is  the  adaptive  detection  of  an  active  sonar  return 
with  robustness  against  doppler,  time  and  frequency  spreading,  and  corrupted  by  reverberation  that 
is  stationary  over  both  data  blocks  and  spatially  distinct  from  the  desired  signal. 

This  structure  may  also  be  applied  to  a  signal  that  has  a  definitive  onset  but  a  long  time  duration. 
This  type  of  signal  suggests  a  sequential  detector  -  this  is  not  the  focus  of  this  memorandum, 
however,  it  may  be  a  future  research  topic. 

3.2  Detector  Derivation 

In  the  Neyman-Pearson  sense  of  maximizing  detection  probability  for  a  fixed  false  alarm  rate, 
substituting  the  observed  data  into  its  likelihood  ratio  yields  the  optimal  detector.  When  the 
statistical  characterization  of  the  data  contains  parameters  with  unknown  values  there  is  often  no 
uniformly  most  powerful  (UMP)  test;  that  is,  a  single  detector  that  provides  optimal  performance 
over  all  possible  values  of  the  unknown  parameters.  The  GLRT  often  provides  good  (although 
suboptimal)  performance  when  there  are  unknown  parameters  and  no  UMP  test  exists. 

In  the  primary-secondary  data  structure  problem,  Kelly  [4]  has  derived  a  GLRT  detector  for  the 
single  frequency  deterministic  signal  case  with  a  single  time  snapshot  PDB.  Robey  et  al.  [5]  and 
Chen  and  Reed  [7]  derived  a  test  by  forming  a  GLR  statistic  for  the  unknown  complex  signal 
amplitude  assuming  that  the  interference  covariance  matrix  is  known.  The  secondary  data  is  then 
used  to  form  the  maximum  likelihood  estimate  (MLE)  of  the  interference  covariance  matrix,  with 
the  estimate  substituted  into  the  above  GLR  statistic.  In  this  memorandum,  a  detector  for  the 
Gaussian  signal  case  is  developed  following  Robey  et  al.  [5]  and  Chen  and  Reed  [7],  noting,  as 
seen  in  [5],  that  the  GLRT  (i.e.,  Kelly  [4])  has  no  assumed  optimality  over  this  method  due  to  the 
nonexistence  of  a  UMP  test. 

The  GLRT  requires  maximization  of  the  likelihood  ratio,  or  equivalently,  the  log-likelihood  ratio, 
with  respect  to  the  unknown  signal  powers,  resulting  in  the  test  statistic 
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_  max 
A = £  i 

m.1  I'm/  j 


log 


f K{Xml\S ml  >Qm) 


(ID 


where  the  summation  indices  m  and  /  represent  frequency  and  time  snapshot  in  the  PDB.  The 
functions  /;(•)  are  the  probability  density  functions  (PDFs)  of  the  primary  data  vectors  under  the 

null  ( H )  and  alternative  ( K )  hypotheses  for  the  Gaussian  signal. 


Substituting  the  N  variate  complex  normal  PDF's  [11]  for  the  Gaussian  signal  model  and 
simplifying  yields  the  form 


mj 

=x 


max 

Sml  . 

max 

•*m/ 


exp[-x"(.v,„,d,„d"  +  Qm)''xm/  +  x",Q;lxm/]| 

■> 

-iog(i+5m/d"QX)'. 

(12) 


Maximizing  this  form  under  the  consu  aint  that  sml  >  0  requires  using  the  signal  strength  estimate 


*m/=i 


1 


l^Q'M 

m  v/ji 


dmQmXm/l‘ 

dmQmd„, 


0 


-1 


Tm,*  1 


^m/  <  1 


(13) 


where  Tml  is  the  normalized  adaptive  beam  output  from  the  m‘h  DFT  bin  and  the  Ith  time  snapshot 
in  the  PDB;  that  is 


Tmi  — 


(14) 


It  is  not  expected  that  there  will  be  signal  energy  in  all  cells  of  the  PDB,  so  it  is  possible  for  the  test 
statistic  A  to  be  corrupted  by  non-signal  cells.  Intuitively,  it  seems  that  if  there  are  a  large  number 
of  non-signal  cells  in  the  PDB,  performance  may  degrade.  A  method,  attributed  to  Wolcin  [1],  for 
dealing  with  this  problem  is  to  form  the  MLE  of  the  cell  signal  power  under  the  constraint. 
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where  the  signal-to-interference  ratio  for  the  ml"1  cell  is  defined  as 

SI*m,  =  Sml^m  Q«i  • 


(15) 


(16) 


This  technique  essendally  fomis  an  estimate  of  where  the  signal  lies  in  the  time-frequency  block  of 
data;  that  is,  if  the  estimated  SIR  is  not  greater  than  t  for  a  given  cell  then  it  is  considered  to  have 
interference  only  and  is  not  incorporated  into  the  test  statistic.  Choosing  t=0  results  in  the 
conventional  MLE  for  the  unknown  signal  amplitude,  since  the  SIR  is  zero  if  the  signal  power  is 
zero.  The  effect  on  the  signal  power  estimate  of  (13)  is 


sm,  =  < 


1 


l"Q"‘d 

m  m 


dwO_1x 

m  ml 

d"Q"‘d 

m  m 


0 


Tm,Z  1+t 
Tml<  1+T 


(17) 


Substituting  this  more  general  estimate  back  into  (12)  yields  the  detector  statistic 


(18) 


where  the  SG  detector  nonlinearity  has  the  following  form: 


gscCn = 


T-l-  log(7") 
0 


r>i  +  r 
T  <  1  +  t' 


(19) 


Substituting  the  complex  Gaussian  PDF  for  the  deterministic  signal  case  into  the  log-likelihood 
ratio  form  of  (1 1)  yields 


max 


*SD 


» ■  ^  ***u/ii  .  \H  \  /  \ 

-jLi„  \  xw^x-.r(xw«A)  (20) 

ml  lUm/J 


By  completing  the  square  in  terms  of  am!,  it  is  seen  that  choosing 
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-  ^  d"Q;,‘x,„; 

dwQ_1d 


will  maximize  the  log-likelihood  ratio  to 


a5D  =  I 


d"Q;‘xmi 


m/  d"Q;X, 

=  8sd{T ml) 


ml 


(21) 


where  the  deterministic  signal  detector  non-linearity  is  actually  the  linear  function 


gso(T)  =  T. 


(22) 


This  is  commonly  known  as  an  energy  detector  and  has  been  explored  for  a  single  frequency, 
single  PDB  time  snapshot  case  by  Robey  [5],  Chen  and  Reed  [7],  and  for  a  time  domain 
implementation  by  Baggenstoss  and  Kay  [12]. 


Substituting  the  MLE  of  the  interference  covariance  matrix  at  the  m'h  DFT  frequency  over  the 
secondary  data  block  at  each  frequency, 

Q.=4iy^y.".  (23) 

J  J*  1 


into  the  adaptive  beam  output  (14)  yields 


dHQ‘ 


km/ 


dH6~‘d 

m  n 


(24) 


3.3  Statistical  Description  of  the  Adaptive  Beam  Output 

The  adaptive  beam  output  (24)  formed  by  using  the  MLE  of  the  interference  covariance  matrix  is 
used  in  both  the  SG  and  SD  detectors.  In  order  to  calculate  the  thresholds  required  to  implement 
the  detectors  and  to  analyze  their  performance,  the  probability  density  function  of  the  adaptive 
beam  output  needs  to  be  determined.  The  approach  used  is  to  first  determine  the  density  function 
of  the  adaptive  beam  output  conditioned  on  the  interference  covariance  matrix  estimate  for  the 
Gaussian  and  deterministic  signals.  The  next  step  is  to  statistically  describe  the  parameters  of  these 
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conditional  distributions  that  depend  on  the  estimated  covariance  matrix.  The  randomness 
introduced  by  estimating  the  interference  covariance  matrix  may  be  described  by  two  independent 
scalar  random  variables,  the  foundation  of  which  comes  from  the  work  of  Reed,  Mallett,  and 
Brennan  [3],  and  is  shown  in  Appendix  A. 

The  adaptive  beam  output,  conditioned  on  the  interference  covariance  matrix  estimate,  is  the 
magnitude  squared  of  a  univariate  complex  Gaussian  random  variable.  As  a  general  case,  consider 
the  N-variate  complex  Gaussian  random  vector 

x  ~  CNN(\i, I), 


and  the  form 


r  =  |wwxf,  (25) 

for  some  constant  vector  w.  The  temporary  variable,  Z  =  wwx,  has  the  density  function  of  a 
univariate  complex  Gaussian  random  variable,  i.e., 

Z~CW(wVww£w). 

The  magnitude  squared  of  Z,  which  is  equivalent  to  (25),  is  the  sum  of  the  squares  of  two 
independent  Gaussian  random  variables  with  the  same  variance.  This  sum  is  proportional  to  a 
non-central  Chi-squared  random  variable  with  two  degrees  of  freedom.  Thus,  when  appropriately 
scaled  by  the  constant  k,  (25)  is  distributed  as 

v^t-xKQ, 

which  has  probability  density  function  [13] 

/r(>’)  =  ^ex  p^-^Y^j/0(Vyf)  y>0, 

where  70(*)  is  a  modified  Bessel  function  of  zero  order.  The  scale  and  non-centrality  parameters 
are,  in  this  case. 
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K  = 


wwIw 


and 


IwV 

C  =  2L77^-. 
wwZw 


(26) 


The  mean  and  power  of  the  adaptive  beam  output,  T,  in  terms  of  the  scale  and  non-centrality 
parameters,  are 


and 


E[T)  =  k{  2  +  C) 

£[72]=  ^(^+8^+8). 


(27) 


The  adaptive  beam  output  (24)  may  be  written  as  in  (25)  by  letting 


w  = 


6"1  d 

^  m  n 


Vd"Q»‘d». 


(28) 


By  substituting  (28)  and  the  mean  and  covariance  from  equation  (8)  or  (9)  into  equations  (26)  and 
(27)  it  is  seen  that  the  chi-squared  scale  and  non-centrality  parameters  and  the  adaptive  beam  output 
mean  and  power  may  be  written  as  functions  of  the  following  two  independent  Gamma  and  Beta 
distributed  random  variables: 


and 


dwO-1d 

**>»  vffi  m 

dwO-1d 

m  ^£#11  m 


Gamma(J  -  N  + 1 ,  j/^) 


Beta(J-N  +  2,W-1). 


(29) 


The  density  functions  for  the  Gamma  and  Beta  random  variables  and  the  derivation  of  these  results 
are  found  in  Appendix  A. 

Utilizing  the  statistical  descriptions  of  the  Gaussian  (8)  and  deterministic  signals  (9),  the  scale, 
non-centrality,  mean  and  power  of  the  respective  adaptive  beam  outputs  are  as  shown  in  table  3-1. 
In  considering  the  results  found  in  table  3-1,  it  is  seen  that  the  Gaussian  signal  case  always  results 
in  a  central  Chi-square  distribution,  which,  with  two  degrees  of  freedom,  is  the  exponential 
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distribution.  The  deterministic  signal  case  results  in  the  exponential  distribution  only  when  the 
SIR  is  zero,  i.e.  when  no  signal  is  present.  The  conditional  mean  of  the  adaptive  beam  output  is 
identical  for  both  signal  types.  The  power  for  the  SD  case  is  always  less  than  or  equal  to  that  of 
the  SG  case  with  equality  only  when  the  SIR  is  zero.  The  only  dependence  of  the  scale  and  non¬ 
centrality  parameters,  that  describe  the  conditional  distribution  of  the  adaptive  beam  output,  on  the 
true  interference  covariance  matrix,  is  through  the  SIR  as  defined  in  (16).  When  no  signal  is 
present,  the  SIR  is  zero  and  the  adaptive  beam  output  distribution  is  free  from  dependence  on  the 
interference  covariance  matrix.  Thus,  any  detector  based  on  the  adaptive  beam  output  will  have  a 
constant  false  alarm  rate. 


Table  3-1  Scale,  non-centrality,  mean  and  power  for  the  adaptive  beam  output 
from  the  ml,h  cell  in  the  PDB  for  Gaussian  and  deterministic  signals. 


Signal 

Scale  k 

Non-cenU'ality  £ 

Mean  E[Tml\ 

Power  £[7y 

SG 

2  pmrm 

0 

pmrm 

■mm 

SD 

1 

2 Pmrm 

2  P„,SIR„,i 

mSSmu 

It  is  useful  to  consider  how  the  Gamma  and  Beta  random  variables  change  with  the  dimension  of 
the  data,  N,  and  the  size  of  the  secondary  data  block,  J .  When  N  -\  the  beamspace  matrix 
preprocessor  is  performing  conventional  beamforming  and  it  is  seen  that  the  Beta  random  variable 
becomes  degenerate  at  one;  that  is,  p  =  1  with  probability  one.  The  secondary  data  is  thus  used  to 
estimate  the  unknown  conventional  beam  output  power  and  then  used  to  normalize  the  primary  data 
as  in  cell  averaging  CFAR  systems.  This  is  a  realistic,  implementable,  non-adaptive  detector.  If, 
additionally,  the  size  of  the  SDB  goes  to  infinity,  which  is  not  a  practical  assumption,  both  the 
Gamma  and  Beta  random  variables  become  degenerate  at  one.  This  situation  represents  infinite 
time  averaging  of  the  unknown  interference  covariance  matrix  which  is  equivalent  to  perfect 
knowledge  about  the  interference  structure.  This  case  is  known  as  the  ideal  conventional  detector 
because  of  the  assumed  knowledge  of  the  interference  statistics  at  the  beam  output  level.  These 
results  are  summarized  in  table  3-2. 
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Table  3-2  Gamma  and  Beta  random  variables  under  certain  N  and  J . 


N 

J 

Gamma  RV 

Beta  RV 

N  >  1 

N  <  J  <  °° 

r,nl 

P  ml 

N  =  1 

N  <  J  <  °° 

r„,l 

Pm,  =  1 

N  =  1 

J  OO 

r,nl  1 

Pm!  1 

4  DETECTOR  ANALYSIS 

In  section  3,  detectors  were  derived  based  on  Gaussian  and  deterministic  signal  assumptions.  The 
resulting  detectors  were  expressed  as  the  summation,  over  all  DFT  frequency  -  time  snapshot  cells 
in  the  PDB,  of  a  univariate  statistic,  the  adaptive  beam  output,  evaluated  at  a  detector  non-linearity. 
The  detection  statistic  is 

A  =  (30) 

ml 


where  £,  (•)  represents  the  detector  non-linearity.  When  the  number  of  cells  in  the  PDB  is  large 
enough,  a  central  limit  theorem  (CLT)  based  approximation  to  a  Gaussian  random  variable  may  be 
used  to  analyze  the  performance  of  these  detectors.  This  requires  knowledge  of  the  mean  and 
variance  of  the  test  statistic.  A,  under  the  no-signal  hypothesis  ( H )  and  the  signal-present 
hypothesis  ( K ).  The  time-frequency  characteristics  of  the  signal  have  not  been  restricted  in 
deriving  the  detectors.  But,  in  order  to  make  analysis  tractable,  certain  restrictions  need  to  be 
placed  on  the  strength  and  time-frequency  location  of  the  signal  in  the  PDB. 

The  signal  restrictions  and  the  mean  and  variance  of  the  detector  test  statistics  are  respectively 
described  in  sections  4.1  and  4.2.  The  detector  performance  analysis  using  the  CLT 
approximation  is  described  in  section  4.3. 

4.1  Signal  Time-Frequency  Model 

Without  some  simplification  of  the  signal  model,  the  analysis  of  the  detectors,  in  particular  the 
variance  computation,  is  difficult  due  to  the  correlation  between  adaptive  beam  outputs,  at  the  same 
DFT  frequency,  introduced  by  the  estimated  interference  covariance  matrix.  A  model  that  provides 
some  freedom  in  the  placement  of  signal  energy  in  the  cells  of  the  PDB  is  desired.  First,  it  is 
assumed  that  when  the  signal  is  present  in  a  time-frequency  cell,  it  will  have  a  constant  SIR.  Next, 
the  number  of  time  snapshot  cells  that  have  signal  present,  for  a  particular  DFT  frequency,  is  either 
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zero  or  Ls.  The  number  of  DFT  frequencies  that  contain  Ls  signal  cells  will  be  Ms.  Thus,  the 
total  number  of  signal  cells  is  MtLs  and  we  define  the  fraction  of  signal  cells  to  total  cells  to  be 


7  = 


AU, 

ML  ' 


(31) 


This  is  graphically  described  in  figure  4-1  with  an  example  where  Ls  =  3,  Ms=  9,  L  =  7 ,  and 
M  =  14.  The  positioning  of  the  DFT  frequencies  that  contain  signal  and  the  signal  cells  at  those 
DFT  frequencies  is  not  restricted. 


M 


Primary  Data  dpt  Bin 

Block  (PDB)  Outputs 


1 


Figure  4-1  Graphical  representation  of  signal  time-frequency  location. 
Signal  is  present  in  3  of  7  time  snapshot  cells  (Ls=  3,  L  =  7)  for  9  of  14 
DFT  frequencies  ( Ms  =9,  M  =  14). 


1  L 

Time  Snapshots  From 
Non-overlapped  DFTs 


4.2  Test  Statistic  Mean  and  Variance 

The  SG  and  SD  detector  functions  were  derived  in  section  3.  In  this  section  the  performance  of  the 
SG  detector  for  the  Gaussian  signal  and  the  SD  detector  for  both  the  Gaussian  and  deterministic 
signal,  as  depicted  in  figure  4-2,  are  analyzed. 


Case 


Deterministic - 

V/aac 

- 1 - 

8so(T) 

Signal 

Detector  Non-linearity 

Gaussian  ■*=— 

- 3 - 

- ►  SscV) 

Figure  4-2  Signal  and  detector  descriptions  for  analysis  cases. 
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The  block  level  mean  and  variance  (i.e.,  mean  and  variance  of  (30))  under  the  signal  present 
hypothesis  ( K )  may  be  found  in  terms  of  the  previously  described  signal  time-frequency 
parameters  and  the  mean  and  variance  of  the  detector  non-linearity,  #,(•),  operating  on  a  single 

adaptive  beam  output,  T,  with  no  signal  and  with  a  signal  with  a  specified  SIR.  These 
expressions  are  given  by 

£*M  =  M.{(  1  -  y)£i[ft(7-)]+  r  £„„[«,  m]} 
and 

VarK[X]  =  Mi(l  -  y)Var0[gj(T)}  +  MLy  Varm[g,(T)\ 

+Mi[i  - 1  +  y(L, -2L+ 1)](«0.0)  -  £„[s,(7-)]!) 

+2MLy(L  -  L'ptO.SIR)  -  £0[s,(7-)]£s„[*,(r)])  <32> 

+M£r(£,  - 1  )(h(,sm,sm  -  £*,,[«, mf) , 

where  the  subscript  on  the  expectation  terms  on  the  light  hand  sides  of  the  equations  denote  the 
SIR  for  the  cell  creating  the  adaptive  beam  output,  T .  Tbe  term  on  tbe  second  Woe  to  \We  Nurvanee 
is  due  to  the  correlation  between  the  adaptive  beam  outputs  at  the  same  frequency.  The 
intermediate  function,  h(a,b),  is  defined  as 

A(a.i)  =  £[&(T,)ft(n)].  (33) 

where  Ta  and  Tb  are  adaptive  beam  outputs  formed  from  the  same  interference  covariance  matrix 
estimate  and  two  independent  primary  data  vectors  with  signal  to  interference  ratio  equal  to  a  and 
b  respectively.  The  mean  and  variance  under  the  no  signal  hypothesis  ( H )  may  be  found  by 
setting  y  =  0  and  Ls  =  0  in  (32),  i.e., 

£K[A]  =  M££,,[si(r)] 

Var„[X]  =  ML  Var„[S,(T)]  +  ML(L - 1)(/>(0,0)  -  £o[s,(D]!)  (34) 


Thus  we  only  need  to  evaluate  the  cell  level  mean  and  variance  (i.e.,  the  mean  and  variance  of  the 
detector  non-linearity)  and  the  function  h(a,b)  in  order  to  determine  the  detector  performance. 
Closed  forms  exist  for  the  SD  detector  for  both  signal  types  and  are  shown  in  table  4-1.  As 
expected,  the  only  difference  between  the  two  signal  types  is  that  the  variance  of  the  adaptive  beam 
output  for  the  deterministic  signal  is  less  than  that  for  the  Gaussian  signal. 
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The  r  and  p  random  variables  in  table  4-1  are,  respectively,  Gamma  and  Beta  distributed,  as 
discussed  in  section  3.3,  and  have  closed  form  inverse  moments  when  J  >  N  + 1.  The  required 
inverse  moments  are 


£[r-']  = 


J 

J-N' 

J 

J-N  +  l' 


j 1 

(j-N)u-N-iy 

JU- 1) 

(J-N)(J-N  + 1) 


(35) 


Table  4-1  Cell  level  mean,  variance  and  h(a,b)  function  for  the  deterministic 
signal  detector  operating  on  deterministic  and  Gaussian  signals. 


Deterministic  Signal 

Detector 

Parameter 

Deterministic  Signal 

Gaussian  Signal 

£[r-‘]{s«+£[p-‘]} 

4r-'}{siR+4p-‘]} 

V°rm[T] 

{5/£  +  £[p"1  ]}3  {2£tr-1  ]  -  ]J| 

{sW  +  £[p-']}J{2£[r-J]-£[r-‘]!} 

+2£[r-1]{£[p-,]-£[p-']I-iS«!} 

+2£[r"3]{£[p'1]-  £[p"]3} 

hSD(a,b) 

E[r~2  }{ab  +  (a  +  b)E[p~]  ]  +  £[p-2  ]} 

E[r-2  ]{ab  +  (a  +  b)E[p^\ +  £[p'2]} 

Due  to  the  non-linear  form  of  the  Gaussian  signal  detector,  the  analysis  is  considerably  more 
difficult  and  is  done  only  for  the  Gaussian  signal  case.  By  utilizing  the  conditional  distribution  of 
the  adaptive  beam  output,  reducing  the  expectation  over  the  Gamma  random  variable,  and  using  the 


transformation  u  - 


1  +  T 


the  ilh  moment  of  the  SG  detector  non-linearity  may  be  expressed  as 


~  Epu 


yK’(i+Ps/R)[,,^r 


(36) 


where  p  is  Beta  distributed  as  previously  described  and  u  is  a  uniform  random  variable 
(independent  of  p )  with  density  function 

f(u)  =  1  0<«<1.  (37) 

The  function  h(a,b)  may  similarly  be  expressed  as 
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hsc{a,b)  =  £pilv 


(J-N+l)(J-N+2)pi{\  +  r)2gM(lii)^(lii) 

jv-vHupaKnpb\^{^£^f' 


(38) 


where  both  u  and  v  are  independent  uniform  random  variables  with  the  density  function  described 
in  (37).  In  reducing  the  expectation  over  the  Gamma  random  variable,  the  following  two 
relationships  are  used: 


and 


Er[rerB]  = 


J-N  + 1 
J(l-e/j)J~N+2 


f  r,VSi -  V-N  +  W-N  +  2) 

rl  J  j\\-%y-N+3 


(39) 


The  evaluation  of  these  functions  would  requite  two-  and  three-dimensional  integrals  for  each  SIR. 
Therefore,  Monte-Carlo  integration  is  used  with  the  uniform  and  Beta  random  variables  generated 
only  once.  The  concept  of  Monte-Carlo  integration  is  discussed  in  Appendix  B. 

4.3  Detection  Performance  via  CLT 

In  order  to  implement  and  evaluate  the  detector,  a  threshold  must  be  chosen  based  on  the  desired 
false  alarm  probability.  The  false  alarm  probability  is  the  probability  of  the  test  statistic  crossing 
some  threshold,  l,  when  no  signal  is  present;  that  is, 

Pra=Pr{A>/|//}.  (40) 


Since  the  test  statistic  is  approximated  by  a  Gaussian  distributed  random  variable  with  mean  and 
variance  as  described  in  (34)  when  no  signal  is  present,  the  threshold  that  provides  the  specified 
false  alarm  performance  may  be  written  as 

/  =  £„[A]-<I>-'(Pn)V^W’  (41> 

where  the  function  $(•)  is  the  standard  Gaussian  cumulative  distribution  function,  defined  by 

1  ii 

<J>(x)  =  Pr{Z  <x}=  e  2  dz,  (42) 
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and  the  superscript  represents  a  functional  inverse. 

The  probability  of  detecting  a  signal  is  the  probability  that  the  test  statistic  crosses  a  threshold  (set 
by  the  desired  false  alarm  rate)  when  the  signal  is  present,  i.e., 

PD  ~  Pr{A  t  l\K].  (43) 


Using  the  CLT  normal  approximation  once  again,  the  probability  of  detection  is 


Po 


■JVarA/i] 


X 

) 


(44) 


The  derivation  of  these  forms  for  the  threshold  and  probability  of  detection  is  straightforward  and 
is  found  in  Nuttall  [2]. 


The  detection  probability  is  dependent  on  the  false  alarm  probability,  the  SIR  found  in  each  signal 
cell,  the  total  number  of  cells  containing  signal,  the  signal  location  parameters  previously 
discussed,  the  adaptive  dimension,  and  the  fixed  sizes  of  the  PDB  and  SDB.  Here,  the  SIR 
required  to  achieve  desired  false  alarm  and  detection  probabilities  is  considered.  Clearly,  the  SIR 
enters  into  the  detection  probability  through  the  mean  and  variance  of  (32),  however,  not  in  a 
simple  fashion.  In  general,  the  SIR  required  to  achieve  specified  detection  and  false  alarm 
probabilities  may  be  found  numerically  by  evaluating  (44)  until  the  desired  detection  probability  is 
found.  The  SG  detector  requires  this  approach,  however,  the  closed  form  solution  to  the  mean  and 
variance  for  the  SD  detector  provides  an  exact  solution.  As  seen  in  table  4-1,  the  cell  level  mean, 
variance,  and  h(a,b)  functions  for  the  SD  detector  are  either  first  or  second  order  polynomials  in 
SIR.  The  block  level  mean  and  variance  are  linear  combinations  of  these  terms  and  may  be 
written,  respectively,  as  a  monomial  and  a  quadratic  in  SIR  as  follows: 

Ek[X]  =  n^SIR  +  rrif,, 

where 

m,  =  MLy£[r-1] 

m1)  =  MLE[r-l]E[p-1]' 
and 

VarK[X]  =  a2SIR 2  +  axSIR  +  a0, 

where 
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a,  =  MLy[i/s£[r-,]  +  i,,(£[r-’]-  £[r-'f )} 

ax  =  2MLy£[p-']|(L+  l)£f;--2]- 1£[/-‘]2} 

a0  =  ML2{E{r-2]E[p-2]  -  £[r-1]2  £[p_1]2}  +  ML£[r'2]£[p"2] 

(45) 


The  Boolean  variable. 


Us  = 


SG 

SD' 


represents  the  difference  between  the  Gaussian  and  deterministic  signals. 

For  the  SD  detector,  we  can  express  the  equality  (44)  in  the  form 

<r'(p0)V^WA]  =  ek[X]  - 1.  (46) 

Following  Nuttall  [2],  we  solve  this  by  squaring  both  sides,  which  yields  a  quadratic  form  in  SIR, 
i.e., 

b2SIR2  +  bxSIR  +  b0  =  0 , 

where 

b2=m?-a2[®-1(PD)]2 

bx  =  2m,  (m^  -  lSD )  -  ax  [O-1  (PD  )]2 ,  (47) 

^o  =  K-4d)2-«o[<i>'1(^)]2 

where  the  threshold  for  the  SD  detector  may  be  written  as 

4d  =  ^-V^’cI>"1(pfJ-  (48) 


The  quadratic  formula  yields  the  solution  to  (47),  however,  only  one  root  will  satisfy  the 
transcendental  equation  (46).  Notice  that  the  above  coefficients  depend  on  the  square  of  4>_1(P0), 
thus,  removing  dependence  on  the  sign.  As  in  [2],  the  sign  of  0-1(P0)  is  used  in  the  quadratic 
formula  to  yield  the  correct  root,  that  is. 
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The  term  in  the  result  may  be  explained  by  allowing  PD  to  increase  slightly  from  one 

half,  which  causes  d>-'(P0)  to  increase  slightly  from  zero,  thus  at  PD  =  0.5,  the  roots  of  the 
quadratic  are  the  same.  An  increased  PD  requires  more  SIR ,  thus  we  want  the  larger  root  which  is 
achieved  in  (49)  when  b2  >  0,  which  is  true  in  the  neighborhood  of  PD  =  0.5. 

As  is  seen  in  Nuttall  [2],  this  method  may  also  be  applied  to  solving  for  the  signal  coverage 
parameter,  y,  required  to  achieve  desired  detection  and  false  alarm  probabilities  using  the  ideal 
conventional  detector  for  a  fixed  SIR.  This  is  done  to  facilitate  plotting  the  required  SIR  against 
the  signal  coverage  parameter.  The  SG  detector  requires  substantial  computation  to  evaluate  the 
required  SIR  for  each  value  of  y.  This  method,  conversely,  computes  the  cell  level  mean  and 
variance  and  finds  the  value  of  y  that  yields  the  desired  performance,  only  requiring  the  numerical 
integration  of  two  equations  for  each  SIR. 

The  ideal  conventional  detector  represents  perfect  knowledge  of  the  conventional  beam  output 
power  which  causes  the  h(a,b)  function  (33)  to  become 

h Ideal  (a'b)  =£.U(j;)te[*(r4  (»» 

because  there  is  no  correlation  between  adaptive  beam  outputs  at  the  same  DFT  frequency  (the  term 
due  to  correlation  in  the  block  level  variance  (32)  thus  goes  to  zero).  Substituting  the  ideal 
conventional  detector  mean  and  variance  into  the  transcendental  equation  (46)  yields  the  quadratic 
form  in  the  signal  coverage  parameter,  y,  i.e., 

c2Y2+cj  +  c0  =0, 

where 

c2  =  MZ.{£OT[«,(D]-£0[i,(D]}2 

c,  =  2 jMLVar^m]  <f'(PrA  ){£,„[*,  (T)]  -  E^r)]} 

- 1®''  (^0)]2{ Var„  [ft  (T)]  -  Var0  [ft  (T)]} 

=  V0r,[ft(D]{[0-(Pra)]! 

As  before,  the  solution  is  the  quadratic  formula  with  a  sign  correction,  that  is, 
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— c,  +sign(*l(pD))^cf  -4c2c0 


2  c, 


(52) 


5  DETECTOR  PERFORMANCE 

In  previous  sections,  Gaussian  and  deterministic  signal  detectors  were  developed  and  a  method  of 
determining  the  SIR  required  to  achieve  specified  detection  and  false  alarm  probabilities  was 
provided.  In  section  5.1,  the  SIR  is  factored  into  an  ideal  case  with  covariance  matrix  estimation 
losses  and  intra-block  correlation  losses  in  order  to  determine  what  has  the  most  effect  on  the 
required  SIR.  It  is  seen  that  the  detector  gain  due  to  adaptive  processing  is  related  to  the  AGI.  In 
section  5.2,  the  AGI  for  Owsley's  standard  test  case  [14]  is  shown,  as  well  as  the  average  loss 
due  to  beam  space  matrix  preprocessing.  In  section  5.3,  the  ideal  conventional  detector  SIR  is 
examined  against  the  fraction  of  cells  containing  signal,  y ,  for  the  SD  detector  with  Gaussian  and 
deterministic  signals  and  for  the  SG  detector  with  various  SIR  thresholds  operating  on  a  Gaussian 
signal.  Estimation  and  correlation  losses  are  then  considered  for  various  size  secondary  and 
primary  data  blocks  and  adaptive  dimension.  From  these  plots,  the  SIR  required  to  achieve 
desired  detection  and  false  alarm  probabilities  can  be  determined. 

5.1  SIR  Decomposition 

It  is  known  that  estimating  the  unknown  interference  covariance  matrix  causes  performance 
degradation  compared  to  the  perfect  knowledge  or  clairvoyant  detector.  These  losses  may  be 
categorized  as  an  estimation  loss  and  as  a  correlation  loss.  The  correlation  loss  results  from  using 
the  same  estimated  covariance  matrix  for  each  time  snapshot  in  the  PDB  at  each  DFT  frequency. 
Analytically,  this  is  represented  by  factoring  the  SIR  into  the  amount  of  SIR  required  for  an  ideal 
conventional  detector,  SIRldcal ,  with  a  loss  due  to  estimation  of  the  interference  covariance  matrix, 
Elou’  and  a  loss  due  to  the  inha-block  correlations,  ;  that  is, 

SIE-Rqd  =  SIR^E^C^,  (53) 

where  the  estimation  and  correlation  losses  are  greater  than  or  equal  to  one.  The  ideal  conventional 
detector  is  a  detector  operating  at  the  output  of  a  conventional  beamfomaer  when  the  beam  output 
power  is  known. 

The  required  cell  level  signal  power,  sRqtl ,  may  be  written,  utilizing  the  definition  of  the  SIR,  as 
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SRi/J  ~  S Ideal 


f  ^Ijoss^Loss  \ 

AG1(N)  )' 


(54) 


where  s,dcal  is  the  required  cell  level  signal  power  for  the  ideal  conventional  detector  for  a  given 
interference  structure  and  AGI(N)  represents  the  AGI  achieved  using  the  P  by  N  preprocessing 
matrix,  Ame.  Thus, 


AG/(yv)=(d"eQl;led„ie)(d»QmaM) 

=  d"A„ie(AfflQmAmfl )  1  A"ad,„  (d"Qmdm ) . 


(55) 


The  loss  due  to  estimation  of  the  interference  statistics  is  defined  as  the  ratio  of  the  required  SIR 
for  a  detector  where  no  correlation  loss  occurs  (L  =  1,  with  M  fixed)  to  that  of  the  ideal 
conventional  detector;  that  is. 


E^J'N) 


SIRj^J,  N,L=l) 


SIR 


Ideal 


(56) 


The  correlation  loss  is  defined  as  the  ratio  of  the  worst  case  required  SIR  (in  terms  of  Ls )  for  the 
given  detector  parameters  to  that  of  a  detector  where  no  correlation  loss  occurs;  that  is. 


max 

r,i  SIRRqd(J,N,L,Lt) 

C.(J,  N,  L)  =  - . 

SIRRqd(J,N,L  =  1) 


(57) 


The  gain  in  adaptive  processing  over  conventional  may  be  expressed  as  the  ratio  of  their  required 
signal  powers,  i.e., 


Gain  = 


1) 


=  AGI(N) 


E,..U,  «  =  D) 

l  Eu,,V,N)  ) 

l  C^J'N'L)  ) 

(58) 


noting  that  AGI(N  =  1)  =  1 .  To  improve  over  the  implementable  conventional  detector,  the  AGI 
need  only  be  greater  than  the  change  in  the  estimation  and  correlation  losses. 
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5.2  Typical  AGI 

As  described  in  (55)  and  applied  in  (58),  the  AGI  represents  the  potential  improvement  when 
adaptive  processing  is  utilized  and  is  a  function  of  the  interference  covariance  matrix  and  the  matrix 
preprocessing.  In  order  to  provide  some  measure  of  what  typical  AGI  values  are,  Owsley's 
standard  test  case  (STC)  [14]  is  analyzed  for  beam  space  preprocessing  on  a  half  wavelength 
linearly  spaced  64  sensor  array.  This  test  case  consists  of  five  plane  wave  signals  amidst  two 
angularly  extended  signals.  In  figure  5- 1 ,  the  expected  conventional  and  adaptive  beam  responses 
to  Owsley's  STC,  Scenario  1,  are  shown.  Any  signal  arriving  from  a  direction  other  than  the  look 
direction  is  considered  to  be  an  interference. 

In  figure  5-2,  the  AGI  achieved  using  element  space  adaptive  processing  is  shown  as  a  function  of 
the  direction  of  arrival  (DOA)  of  the  signal.  The  largest  AGI  values  are  located  close  in  arrival 
angle  to  the  strong  plane  wave  interferences  and  are  seen  to  be  greater  than  10  dB.  In  figure  5-3, 
the  average,  over  arrival  angles,  of  the  loss  in  AGI  incurred  by  the  beamspace  preprocessing  is 
plotted  against  the  adaptive  dimension  of  the  preprocessor  (number  of  beams  used  in  the 
preprocessing  matrix).  This  curve  provides  an  idea  of  how  much  AGI  is  lost  due  to  beamspace 
preprocessing.  In  this  example,  five  beams  provides  an  average  of  0.5  dB  of  loss. 


sin(DOA) 

Figure  5-1  Expected  conventional  and  adaptive  beam  responses  to  Owsley's 
STC,  Scenario  1 ,  for  a  half  wavelength,  linearly  spaced,  64  sensor  array.  The 
marks  on  the  abscissa  indicate  signal  placement. 
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Figure  5-2  Array  gain  improvement  ( AGI )  for  Owsley's  STC,  Scenario  1,  for  element  space 
processing  of  a  half  wavelength,  linearly  spaced,  64  sensor  array. 


Average  AGI  Loss  due  to  Beamspace  Preprocessing 


Figure  5-3  Average  loss  in  AGI  compared  to  element  space  adaptive  processing 
for  Owsley's  STC,  Scenario  1,  due  to  beamspace  preprocessing  as  a  function  of  the 
number  of  beams  used  in  the  preprocessor,  for  a  half  wavelength,  linearly  spaced, 
64  sensor  array. 
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5.3  Performance  Results 

The  method  of  analysis  for  the  SD  and  SG  detectors  used  here  will  first  consider  the  SIR  required 
to  achieve  specified  detection  and  false  alarm  probabilities  for  the  ideal  conventional  detector 
against  the  parameter  y ,  the  fraction  of  cells  in  the  PDB  that  contain  signal.  The  estimation  loss  is 
then  determined  for  a  narrowband  (y  near  zero)  and  a  broadband  (y  near  one)  signal,  with 
changing  SDB  sizes  and  adaptive  dimensions.  Finally,  the  correlation  loss  is  determined  as  a 
function  of  the  PDB  size  for  a  fixed  SDB  size.  Throughout  this  analysis,  the  total  number  of  cells 
in  the  PDB,  ML,  is  held  constant.  The  SD  detector  is  considered  for  both  the  deterministic  and 
Gaussian  signals.  The  SG  detector  is  considered  only  for  the  Gaussian  signal.  These  signal- 
detector  combinations  will  be  denoted  by:  SD-SD,  SG-SD,  and  SG-SG,  with  the  first 
variable  representing  the  signal  type  and  the  second  the  detector  type. 

The  SIR  required  to  achieve  a  PD  =  0.9  and  PFA  =  10"\  for  an  ideal  conventional  detector  with 
ML-  1024,  is  shown  in  figure  5-4  for  the  SD-SD  and  SG-SD  combinations  and  for  the 
SG  -  SG  combination  with  r  =  0, 2,  and  4.  This  plot  corresponds  to  those  found  in  [2]  for  the 
Gaussian  signals. 

For  small  y  signals,  i.e.  signals  that  are  confined  to  very  few  cells  in  the  PDB,  the  SG  detector 
performs  better  with  increasing  threshold,  T.  Raising  the  threshold  too  high,  thus  decreasing  the 
number  of  non-zero  values  in  (18),  violates  the  assumptions  required  to  use  the  CLT  in  the 
analysis  for  small  y  signals.  At  values  of  y  close  to  unity,  i.e.  signals  that  almost  completely 
cover  the  PDB,  the  SD  outperforms  all  of  the  SG  detectors  for  the  Gaussian  signals.  This 
characteristic  is  expected  and  has  been  shown  before  (see  [2]).  An  intuitive  rationale  is  that  the 
deterministic  signal  should  be  easier  to  detect  than  the  Gaussian  signal  because  the  signal  structure 
is  known  exactly.  This  is  demonstrated  in  figure  5-4  for  the  SD  detector  with  the  most  degradation 
for  the  Gaussian  signal  when  y  is  small. 

When  the  interference  statistics  are  unknown  and  need  to  be  estimated,  the  increase  in  the  required 
SIR  is  quantified  by  the  metric.  For  a  fixed  number  of  PDB  cells  ( ML  =  1024),  this 

parameter  is  considered  against  increasing  values  of  SDB  size  (7),  for  N  =  1,5, 10,  and  15  for 
both  small  y  and  large  y  signals.  Figures  5-5  and  5-6,  respectively,  show  the  estimation  loss  for 
the  y  =  0.1  and  y  =  0.9  signals.  The  ideal  conventional  detector  required  SIR  is  noted  in  figure  5- 
4  for  both  of  these  signal  cases.  Recall  that  the  number  of  time  snapshots  in  the  PDB  is  one 
( L  =  1)  so  there  is  no  loss  due  to  con-elation. 
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If  -  Fraction  of  cells  containing  signal 

Figure  5-4  SIRldeal  vs.  y  for  ML  =  1024,  PD  =  0.9,  and  PFA  =  10"3. 


As  expected,  the  loss  due  to  estimating  the  interference  statistics  decreases  as  the  size  of  the  SDB  is 
increased.  The  estimation  loss  seems  invariant  to  y  for  the  SD  -  SD  and  SG  -  SD  combinations, 
with  the  loss  nearly  indistinguishable  for  the  two  signal  types.  The  SG  detector  ( T  =  0)  shows 
increased  loss  for  small  y ,  with  a  higher  sensitivity  for  larger  N .  Note  also  that  the  losses  for 
each  detector  are  approximately  the  same  for  the  large  y  signal. 

In  order  to  consider  the  loss  due  to  correlation  of  the  normalized  beam  output  statistics, 
,N ,L) ,  we  choose  a  particular  size  for  the  SDB  (  /  =  20).  In  most  cases,  the  SDB  size  will 

be  dictated  by  the  stationarity  of  the  interference,  thus  this  is  a  reasonable  specification.  Figures 
5-7  and  5-8  show  the  loss  in  required  SIR  due  to  correlation  against  the  PDB  size  (L)  for 
N  =  1, 5, 10,  and  15  and  respectively  for  the  y  =  0.1  and  y  =  0.9  signals,  respectively. 

The  correlation  losses  for  all  the  signal-detector  combinations  increase  as  the  number  of  time 
snapshots  in  the  PDB  increases.  As  with  the  estimation  loss,  the  loss  for  the  SD-SD  and 
SG  D  combinations  are  indistinguishable.  The  correlation  loss  is  also  seen  to  be  larger  for  the 
small  y  signal  for  all  the  signal-detector  combinations  and  slightly  sensitive  to  larger  N. 
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J  -  Number  of  time  snapshots  in  Secondary  Data  Block 

Figure  5-5  EUss(J,N)  vs.  J  for  y  =  0.1 ,  ML  =  1024,  PD  =  0.9,  and  PFA  =  10'3. 


J  -  Number  of  time  snapshots  in  Secondary  Data  Block 

Figure  5-6  £^(7, N)  vs.  J  for  y  =  0.9,  ML  =  1024,  PD  =0.9,  and  PFA  =  10"3. 
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PFA  =  io-3. 


and  PFA  =  10"3. 
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6  CONCLUSIONS 

This  memorandum  considers  the  utilization  of  multiple  sensor  array  data  in  the  adaptive  detection 
of  sonar  and  radar  signals.  The  signal  processing  prior  to  the  detector,  including  a  dimension 
reducing  matrix  preprocessor  and  statistical  models  commonly  used  to  describe  the  data,  were 
introduced  as  a  preface  to  the  adaptive  detector  development.  Detectors  were  then  developed  for 
the  Gaussian  and  deterministic  signal  cases,  where  signal  energy  is  restricted  to  a  block  in  time  and 
frequency,  the  primary  data  block  (PDB).  The  Gaussian  signal  detector  was  derived  allowing  the 
signal-to-interference  ratio  (SIR)  that  is  observed  at  each  time  snapshot-DFT  frequency  cell  in  the 
PDB  to  be  thresholded.  This  imparts  a  robust  character  to  the  detector  against  diverse  time- 
frequency  characteristics  of  a  received  signal.  A  second  block,  disjoint  in  time  from  the  block 
containing  the  signal,  is  used  to  estimate  the  interference  statistics  which  are  common  to  both 
blocks.  Applications  of  this  data  structure  include  active  sonar  returns  and  limited  time  duration 
signals  that  are  corrupted  by  interference  that  is  stationary  over  both  data  blocks. 

The  detection  and  false  alarm  probabilities  of  the  detectors  were  then  determined  using  a  central 
limit  theorem  based  approximation  to  a  normal  random  variable.  Closed  form  solutions  were 
found  for  the  deterministic  signal  detector  operating  on  both  the  deterministic  and  Gaussian 
signals.  The  Gaussian  signal  detector  required  the  evaluation  of  multi-dimensional  integrals  which 
was  accomplished  by  a  Monte-Carlo  method. 

The  performance  of  the  detectors  was  then  described  by  the  cell  level  SIR  required  to  achieve 
specified  detection  and  false  alarm  probabilities.  The  SIR  was  factored  into  the  SIR  required  for 
an  ideal  conventional  detector  affected  by  an  estimation  loss  and  an  intra-block  correlation  loss. 
When  the  matrix  preprocessing  of  the  array  data  provides  the  detector  with  multivariate  data,  an 
adaptive  processing  gain,  quantified  by  the  array  gain  improvement  (AGI),  may  counter  the 
estimation  and  congelation  losses.  The  AGI  increases  with  the  dimension  of  the  preprocessed  data, 
when  beamspace  preprocessing  is  used,  however,  this  improvement  is  opposed  by  an  increase  in 
the  estimation  and  correlation  losses.  The  estimation  loss  decreases  as  the  time  length  of  the  block 
used  to  estimate  the  interference  statistics  increases.  The  correlation  loss  increases  with  the  time 
length  of  the  block  containing  the  signal.  Both  the  estimation  and  correlation  losses  are  seen  to 
increase  with  the  adaptive  dimension  of  the  detectors. 

These  results  lead  to  a  desired  detector  structure  that  has  (i)  a  large  time  block  used  to  estimate  the 
interference  statistics,  (ii)  a  small  time  block  that  contains  the  signal,  (iii)  a  small  adaptive 
dimension  to  reduce  estimation  and  correlation  losses,  and  (iv)  a  large  adaptive  dimension  to 
improve  the  AGI.  Clearly,  (iii)  and  (iv)  are  conflicting  requirements,  however,  the  adaptive 
dimension  may  be  reduced  substantially  using  beam  space  preprocessing  with  a  loss  in  AGI  that  is 
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easily  countered  by  the  reduced  estimation  and  correlation  losses.  The  large  time  desired  for 
estimation  of  the  interference  statistics,  (i),  is  opposed  by  the  limited  stationarity  of  typical  sonar 
interferences.  The  size  of  the  time  block  containing  signal  is  determined  by  the  type  of  signal  being 
detected  and  by  the  size  of  the  DFT  that  decorrelates  the  time  series  data.  Fewer  time  snapshots  in 
the  block  containing  signal  is  desired  which  requires  a  larger  DFT  which  is  good  for  decorrelating 
the  time  series,  however  it  also  means  fewer  time  snapshots  in  the  block  used  to  estimate  the 
interference  statistics,  thus  increasing  the  estimation  loss.  The  decrease  in  correlation  loss, 
however,  is  not  enough  to  offset  the  increase  in  estimation  loss,  thus  the  DFT  size  should  be 
chosen  large  enough  to  decorrelate  the  time  series  data  and  small  enough  to  provide  adequate 
estimation  of  the  interference  statistics. 

This  memorandum  proposes  two  adaptive  detector  structures  based  on  deterministic  and  Gaussian 
signals  that  exist  in  a  limited  region  in  time  and  frequency.  The  detectors  estimate  and  then  employ 
the  unknown  interference  covariance  matrix  which  results  in  a  constant  false  alarm  rate  detector.  A 
method  for  system  level  design  is  provided  by  evaluating  the  effect  of  the  size  of  the  primary 
(signal  and  interference)  and  secondary  (interference  only)  data  blocks  on  the  SIR  required  to 
achieve  specified  performance  for  several  adaptive  dimensions.  The  size  of  the  DFT,  the  type  of 
signal  to  be  detected,  and  the  length  of  the  stationarity  of  the  interference  affect  these  sizes,  and, 
thus,  the  estimation  and  correlation  losses  incurred  on  the  required  SIR.  Balancing  these  losses  is 
the  gain  due  to  adaptive  processing  which  increases  with  the  adaptive  dimension  of  the 
preprocessing.  Coupling  the  estimation  and  correlation  losses  with  the  adaptive  processing  gain, 
as  a  function  of  the  adaptive  dimension,  an  optimal  choice  of  adaptive  dimension  is  possible  for  a 
given  interference  environment. 
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APPENDIX  A  -  PDF  OF  CFAR  AND  SIR  LOSS  RANDOM  VARIABLES 
The  distribution  of  the  adaptive  beam  output  for  the  m"‘  DFT  frequency  and  the  Ith  time  snapshot 
in  the  PDB  (24)  conditioned  on  the  estimate  of  the  interference  covariance  matrix  was  described  in 
terms  of  SIR  loss  and  CFAR  loss  parameters.  The  SIR  loss  parameter,  was  first  derived  and 
shown  to  be  Beta  distributed  by  Reed,  Mallet  and  Brennan  [3];  its  definition  and  distribution  are 


and 


(d,"QXf 

pm  ~  Beta{J  -  N  +  2,  N  - 1), 


/p(P) 


rv+i) 


J-N+l 


r(j-N+2)r{N-iy 


(]-p) 


N-2 


0<p<  1, 


(A-l) 


where  the  Gamma  function  is  defined  as 


T  (z)=  \xl-xe~x(h c. 

jr=0 


For  positive  integer  values  of  z,  the  Gamma  function  is 

T(z)  =  (z-1)!  z  =  1,2,3,... 


In  this  appendix,  it  will  be  shown  that  the  CFAR  loss  parameter  is  a  Gamma  distributed  random 
variable  and  that  it  is  statistically  independent  from  the  SIR  loss  parameter.  As  seen  in  (29),  the 
CFAR  loss  parameter  is  defined  by 


r 


in 


dwO~‘d 

_  ni 

dw6  ‘d  ' 

m  ^ Cm  m 


(A-2) 


It  is  now  convenient  to  introduce  the  complex  Wishart  distribution  of  JQm : 

A  =  /Q.  -CW*(y,Q„).  (A-3) 

Equation  (A-3)  reads:  the  N-by-N  matrix  A  has  a  complex  Wishart  distribution  with  J  degrees 
of  freedom  and  scale  mauix  Qffl .  Letting  v  =  Q^dnl ,  the  CFAR  loss  parameter  becomes 
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r.„  = 


v"v 


J  vwQ*  A^v 

>Cm  ’ 


I  y 

Letting  B  =  Qm* AQmJ  and  u  =  -j  (so  that  uwu  =  1),  then 

Vvwv 


r..  = 


1 


J  uwB_1u 


(A-4) 


(A-5) 


Let  u  =  Pwe,  and  C  =  PBPW,  where  P  is  orthogonal  (PWP  =  i)  and  ej  is  the  first  unit  vector  (i.e. 
ei  =[1"  Of),  giving 


Now  it  is  seen  that  C  is  a  Wishart  distributed  matrix  by  using  the  property  that  if  A  ~  CWN(J,I i) 
and  M  is  L-by- N  with  rank  L,  then  MAM"  ~  This  is  found  for  the  real 

Wishart  matrix  in  Muirhead  [15]  and  is  easily  extended  to  the  complex  Wishart  case.  Thus  writing 
C  as 


’inW 


C  =  PQm'AQjP 


it  is  seen  to  be  complex  Wishart  distributed, 

C-CtV^.PQ^Q.Q-jp") 


In  this  same  fashion,  the  SIR  loss  parameter  may  be  written  as 


Pn, 


Proceeding,  we  first  partition  the  maUix  C  and  D  =  C“‘  as  follows: 


(A-7) 


(A-8) 
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and 

(A-9) 

where  cu  and  dn  are  scalars,  C22  and  D22  are  (N-l)-by-(N-l)  matrices,  and  the  other 
dimensions  follow  accordingly.  Now,  clearly, 

e,  C  e,  =  dn 
and 

ei  C  e,  =  dn  +d12d21. 


D  =  C‘  = 


du  < 

|d21  D22 


Either  from  equating  CD  =  I  or  from  the  Matrix  Theory  Appendix  of  [15], 


and 


d\\—{p\\  C12^-'22C2l) 


-C"‘c 

d  = - — 22_21 


-l 


c  -cHC~lc 

C,,  «-i2v'22i'21 


(A- 10) 


:  dx  iC22c2, 


The  complex  Wishart  matrix  form  of  a  theorem  found  in  Muirhead  ([15]  p.93  Thm.  3.2.10)  is  now 
required.  The  extension  from  the  real  Wishart  matrix  form  is  simply  performed  by  following 
Muirhead's  proof  and  using  the  complex  Wishart  distributional  information  found  in  [11]. 

Theorem  -  Suppose  that  A  ~  CWA,(7,Z),  with  the  partitions 


A  = 

'A„ 

A12 

and  E  = 

S.2 

_A21 

A22. 

*2, 

^22. 

where  Au  and  En  are  M-by-M  matrices.  Now  define 


Aj  i*2  Aj  |  A|2A22A21 

and 

^11-2  =  ^11  ~  ^12^22^21’ 

then 

O')  A u.2~CWm{J-N+M,Lu.2) 

(ii)  Au,2  is  independent  of  A12  and  A22 

Now,  letting  M  - 1 ,  A  =  C,  and  E  =  /  in  the  above  theorem,  it  is  seen  that 
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and 


e,WC-1e1=(C11,r 

ef,Cae1=(l  +  c*<5cJ1Xc11.a)'a. 


(A-ll) 


Substituting  these  into  (A-6)  and  (A-8)  yields 


r  _  fil-2 
"  J 


and 


Pm  = 


1 


1  +  c.^ 

i  i  ^12v*x22^'21 


(A-12) 


Utilizing  the  above  theorem,  it  is  seen  that  rm  and  pm  are  independent  and  that  the  distribution  of 
Cu,2  =  Jrm  is  CW,(7-/V  +  1,1),  which  is  a  Gamma(J-N  + 1,1)  distribution,  thus  rm  is 
distributed  as 


with  PDF 


rm  ~  Gamma(J  -  N  + 1,1/7), 


iJ-N+l 

fr(r)=_~-  r  >  0. 


r(7-yv  +  i) 


(A- 13) 
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APPENDIX  B  -  MONTE-CARLO  INTEGRATION 

It  is  very  difficult  to  evaluate  higher  dimension  integrals  numerically  [16].  To  avoid  this  problem, 
Monte-Carlo  integration  [16]  is  employed  for  evaluating  the  two  and  three  dimensional  integrals  in 
this  report.  This  method  involves  writing  the  integral  as  the  expectation  of  a  function  over  a  set  of 
(in  this  case)  independent  random  variables, 

'¥  =  \g(x)dF(x)  =  Ex[g(x)l  (B-l) 


where  g(\)  is  the  function  to  be  integrated  and  F(x)  is  the  joint  cumulative  distribution  function  of 
the  random  vector  x . 

The  integral  is  then  approximated  by  forming  a  large  number  of  these  random  variables  and 
computing  the  sample  mean  of  the  desired  function, 

@-2) 

M  n=l 

where  the  x„  are  independent  samples  of  the  random  vector  x.  This  is  known  to  converge  in 
probability  to  the  true  value  of  the  integral  as  the  number  of  samples  used  tends  to  infinity. 
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